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Abstract 

The phenomenon of stochastic eternal inflation is studied for a chaotic inflation 
potential in a Bianchi Type I spacetime background. After deriving the appropriate 
stochastic Klein-Gordon equation, we give details on the conditions for eternal inflation. 
It is shown that for eternal inflation to occur, the amount of anisotropy must be small. 
In fact, it is shown that eternal inflation will only take place if the shear anisotropy 
variables take on values within a small region of the interior of the Kasner circle. We 
then calculate the probability of eternal inflation occurring based on techniques from 
stochastic calculus. 


1 Introduction 

The theory of inflation arose out of trying to explain the so-called flatness and horizon prob¬ 
lems in cosmology [T]. However, the vast majority of non-eternal and eternal inflationary 
models assume from the outset that the background spacetime is spatially homogeneous, 
flat, and isotropic, and are therefore considered in a /c = 0 Friedmann-Lemaitre-Robertson- 
Walker (FLRW) spacetime [2]. Stochastic eternal inflation has received a considerable 
amount of attention over the past few years as it is one of main motivations behind the 
multiverse [2]. In this paper, we study the dynamics of stochastic eternal inflation in the 
context of a Bianchi Type I spacetime, which is spatially homogeneous, but anisotropic, 
and hence, is more general than the standard k = 0 FLRW cosmology. It is of considerable 
interest to consider whether stochastic eternal inflation could happen in a more general 
background spacetime as this would allow the possibility of the existence of a multiverse 
to arise from more generic conditions, in addition to solving the so-called “fine-tuning” 
problem in a more general sense. To the best of the authors’ knowledge, such a study has 
not been presented in the scientific literature. 
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Non-eternal inflationary dynamics in the context of Bianchi Type I cosmological models 
have been studied a number of times in the literature. Barrow and Turner [3] showed that 
a universe with moderate anisotropy will undergo inflation and will rapidly be isotropized. 
Rothman and Marsden [3] analyzed in detail inflationary dynamics in Bianchi Type I cos¬ 
mologies. Jensen and Stein-Schabes [5] studied the effects of anisotropic cosmologies on 
inflation in significant detail. Moss and Sahni [6] studied the effect of anisotropies on 
the chaotic inflationary for a general class of Bianchi cosmologies. Rothman and Ellis [7] 
showed that inflation takes place in all Bianchi models and that anisotropy does not sup¬ 
press inflation. Belinsky and Khalatnikov [8] studied homogeneous cosmological models of 
Bianchi I and Friedmann type in the presence of a massive scalar field. They showed that 
the majority of solutions undergo an inflationary stage. Grpn ^ considered inflationary 
Bianchi Type I models with shear, bulk, and nonlinear viscosity. Kitada and Maeda m 
proved a cosmic no-hair theorem for Bianchi models in power-law inflation. Parnovskij m 
analyzed the evolution of several Bianchi Type models in the presence of a scalar field. R 
was shown that certain Bianchi models reduce to the Bianchi Type I model in this case. 
Burd and Barrow [12] considered cosmological models with a scalar field with an exponen¬ 
tial potential. They investigated the inflationary nature of a re -|- 1-dimensional Friedmann 
universe and some 3 -|- 10-dimensional cosmological models as well. Aguirregabiria, Fein- 
stein, and Ibanez m obtained a general exact solution of the Einstein field equations for 
Bianchi type I universes filled with an exponential-potential scalar field and also studied 
their dynamics. Ibanez, van den Hoogen, and Coley m studied whether homogeneous 
cosmological models containing a self-interacting scalar field with an exponential poten¬ 
tial isotropize. Harko m studied the dynamics of a Bianchi type I universe filled with 
a scalar field, and also discussed conditions under which inflation can occur in such mod¬ 
els. Chimento m found the exact general solution for cosmological models arising from 
the interaction of the gravitational field with two scalar fields in both flat FLRW and the 
locally rotationally symmetric Bianchi Type I spacetimes filled with an exponential poten¬ 
tial. Byland and Scialom m studied the Einstein-Klein-Gordon equations for a convex 
positive potential in a Bianchi Type I, a Bianchi Type III, and a Kantowski-Sachs uni¬ 
verse. After analyzing the inherent properties of the system of differential equations, they 
studied the asymptotic behaviors of the solutions and their stability for an exponential 
potential. Chakraborty and Paul [18] studied the cosmic no-hair theorem for anisotropic 
Bianchi models which admit an inflationary solution with a scalar field. It was found that 
the form of the potential does not affect the evolution in the inflationary era. Chimento 
m showed that inflation can be obtained from nonaccelerated scenarios by a symmetry 
transformation in Bianchi Type I spacetimes. Bah and Jain [20] investigated a Bianchi 
Type I inflationary universe in the presence of massless scalar field with a flat potential, 
and discussed the inflationary scenario in detail. Fay and Luminet [21] studied in detail 
different conditions leading to the isotropization of a Bianchi Type I model. Folomeev [22] 
studied the dynamics of a Bianchi Type I model with two interacting scalar helds. Gum- 
rukcuoglu, Contaldi, and Peloso [23] studied inflationary perturbations in a Bianchi Type 
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I model and studied their imprint on the cosmic microwave background. Pitrou, Pereira, 
and Uzan [24] investigated the predictions of an inflationary phase starting from a Bianchi 
Type I universe. Beyer and Escobar [25] considered the dynamics of Bianchi A scalar field 
models which undergo inflation, and in particular, studied the “graceful” exit problem. 

Stochastic eternal inflation has also been studied a number of times in different contexts. 
Mijic [26| analyzed the boundary conditions of the dynamics of inflation as a relaxation 
random process and gave a simple proof for the existence of eternal inffation. Salopek 
and Bond m studied nonlinear effects of the metric and scalar fields in the context of 
stochastic inflation. Linde, Linde, and Mezhlumian [28| considered chaotic inflation in 
theories with effective potentials that behave either as </>"' or as They also performed 
computer simulations of stochastic processes in the inflationary universe. Linde and Linde 
[29] investigated the global structure of an inflationary universe both by analytical methods 
and by computer simulations of stochastic processes in the early universe. Susperregi and 
Mazumdar m considered an exponential inflation potential and showed that this theory 
predicts a uniform distribution for the Planck mass at the end of inflation, for the entire 
ensemble of universes that undergo stochastic inflation. Vanchurin, Vilenkin, and Winitzki 
m investigated methods of inflationary cosmology based on the Fokker-Planck equation of 
stochastic inflation and direct simulation of inflationary spacetime. Winitzki [32] explored 
the fractal geometry of spacetime that results from stochastic eternal inflation. Kunze [33] 
considered chaotic inffation on the brane in the context of stochastic inflation. Winitzki 

[34] described some issues regarding the time-parameterization dependence in stochastic 
descriptions of eternal inflation. Gratton and Turok |35| investigated a simple model of X(j)^ 
inflation with the goal of analyzing the continuous revitalization of the inflationary process 
in some regions. Li and Wang [36] used a stochastic approach to investigate a measure 
for slow-roll eternal inflation. Tomas Galvez Ghersi, Geshnizjani, Piazza, and Shandera 
m used stochastic eternal inflation to analyze the connection between the Einstein equa¬ 
tions and the thermodynamic equations. Qiu and Saridakis [38] used stochastic quantum 
fluctuations through a phenomenological, Langevin analysis studying whether they can 
affect entropic inflation eternality. Harlow, Shenker, Stanford, and Susskind [39] described 
a discrete stochastic model of eternal inflation that shares many of the most important 
features of the corresponding continuum theory. Vanchurin [JO] developed a dynamical 
systems approach to model inflation dynamics. Eeng, Li, and Saridakis m investigated 
conditions under which phantom inflation is prevented from being eternal. Recently, Kohli 
and Haslam [42] showed that the stochastic eternal inflation scenario in the context of a 
spatially flat FLRW cosmology leads to strong singularities, thus eliminating the possibility 
of inflation being eternal in such models. 

In this paper, we consider the effects of adding a stochastic forcing term in the form 
of Gaussian white noise to the right-hand-side of the Klein-Gordon equation as is done in 

[35] and [42] . Following [2], we note that the classical dynamics of the inflaton dictates 
that the inflaton always rolls down its potential. However, quantum fluctuations can also 
drive the inflaton uphill, which causes inflation to last longer in some regions, which, in 
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turn, enlarges the volume of the region. In some regions, the inflaton will remain high 
enough up the potential hill to maintain acceleration. This is a stochastic scenario that is 
typically described as eternal inflation which is one of the main motivating ideas behind the 
multiverse concept. As discussed in [35], the purpose of this stochastic term is to describe 
the scalar field fluctuations. 

Note that, throughout, we assume geometrized units, c = SttG = 1, so that all dynam¬ 
ical quantities have units of powers of length. 


2 The Dynamical Equations 


In this paper, we are considering a Bianchi Type I universe. We will in this section, 
employ the orthonormal frame formalism [33| to the Einstein field equations to obtain a 
dynamical system of equations. Following |44] . we consider a group-invariant orthonormal 
frame {n, eQ,}, where n is the unit normal to the group orbits. Further, since n is tangent 
to a hypersurface-orthogonal congruence of geodesics, all dynamical variables depend only 
on time, t. The Einstein field equations take the form 

H = + (1) 

3 D 

ffafS — ( 2 ) 

p = (3) 


where H is the Hubble parameter, aap is the shear tensor, and ^ and p are the energy 
density and pressure of the matter content in the cosmological model respectively. 

As mentioned in the introduction, our universe model contains a scalar field, which is 
minimally coupled and has the Lagrangian density 


C 


-Vd 




(4) 


which implies that the energy-momentum tensor takes the form 


= V“(/>VV - 


(5) 


Following |2], since we are working within a 3 -|- 1 formalism, such that we assume 
that Va</> is timelike, and is normal to the space like surfaces = const, the energy- 

momentum tensor has the algebraic form of a perfect fluid, with 


At = + 4^(0), P= -y{(!>)■ 


( 6 ) 


The Klein-Gordon equation which is implied from the divergence-free property of Tab 
above then becomes 

4> + 3H(^ + V'icP) = 0. (7) 
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To model the quantum fluctuations in the eternal inflation scenario, we add a stochastic 
noise forcing term r](t) on the right-hand-side of this equation, so that it takes the form 

^ + 3H<j) + V’i^) = ( 8 ) 


Further, we also note that since the shear tensor is diagonal and trace-free, we 
define ^ ^ 

<7+ = - (i722 + 0 - 33 ) , (T_ = ——= ( 1 T 22 — (T 33 ) , (9) 


such that 


2^3 


<7^ = = 3{al + a'i) . 


( 10 ) 


In this paper, we will specifically consider the case of a chaotic inflation potential, 
V{(l)) = and will also make the slow-roll approximation, where V{<i)) > Within 

the slow-roll approximation, the Friedmann equation (l3|), takes the form 


!/(</)) = 3H^ -a^. 

We will also employ expansion-normalized variables, defining, 

^2 






x^.lvcp) 1 


3 H 


such that 


= 1 . 


Under the slow-roll approximation, the Klein-Gordon equation 


becomes 


( 11 ) 

( 12 ) 

(13) 

(14) 


Letting U((/>) = and employing the expansion-normalized variables in Eq. (I12p . 

one obtains 

where we have defined 


X' - {1 + q)X + mix = J\f{T), 


(15) 


W(r) = 


V 


m 


3Vm' 

Note that mo is the expansion-normalized dimensionless inflaton mass, which in our units 
has dimensions of H |45j . 

Further, q is the deceleration parameter which is found from Raychaudhuri’s equation 
(flj) to be 

q = 2Y?-X‘^. (17) 

Finally, the expansion-normalized shear propagation equations take the form 


(16) 


S'± = S±((?- 2 ). 


(18) 
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Note that, in the expansion-normalized variables, differentiation is denoted by a primed 
symbol, and is with respect to a dimensionless time variable r, which is defined as 1441 
dtjdr = IIH. 

The Gaussian noise term M{t) deserves some explanation, which we now provide and 
is based on arguments given in [l2]. We will actually make use of the fact that the noise 
term is formally defined as the time derivative of the Wiener process W{t), that is, 

W'{t) = W(r), W{0) = 0. (19) 

It should be noted that this definition defines the Gaussian noise term as the formal deriva¬ 
tive of the Wiener process. As we show below, the Wiener process is actually nowhere 
differentiable, so one has a set of stochastic differential equations as opposed to ordinary 
differential equations. As noted in [35], M{t) also satisfies the autocorrelation relation 

(W(t)W(t')) = 5{t — t'). (20) 

Further, it is important to realize that W is nowhere differentiable. In fact, for com¬ 
pleteness, following [3^ we now state some properties of W based on the so-called Levy 
characterization: 


1. The path of W is continuous and starts at 0, 

2 . VF is a martingale and [dVF(r)]^ = dr, 

3. The increment of W over time period [s,t] is normally distributed, with mean 0 and 
variance (r — s), 

4. The increments of W over non-overlapping time periods are independent. 


Note that, we will not go into extensive detail about martingales. The interested reader is 
asked to consult m for more details on martingales and their properties. Based on the 
arguments provided in [46|, we will now show that W is not anywhere differentiable, that 
is, not anywhere. Consider a time interval of length At = 1/n starting at r. We will 
define the rate of change over an infinitesimal time interval [r, r -|- Ar] is 




W{t + Ar) - W{t) 
At 


IT(r+i-IF(r) 

1 

n 


w 



W{t) 


( 21 ) 


Therefore, A„ is normally distributed with expectation value 0, variance n, and of course, 
standard deviation, ^/n. It is clear then that Xn has the same probability distribution 
as where Z is the standard normal distribution. To analyze the properties, we 

must see what happens to Xn as Ar —)• 0, in other words as n ^ oo. For any /c > 0, let 
Xn = 'JnZ ^ then 


P[||A„|| >k]=F 


1^11 > 


V^. 
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( 22 ) 








Clearly, as n —oo, kj^fn —0, so we have that 

^ P[||Z|| > 0] = 1. (23) 

The point is that we can choose k to be arbitrarily large, so that the rate of change at 
time r is not finite, and therefore, W is not differentiable at r. Since r is arbitrary, one 
concludes that W is nowhere differentiable. 

In summary, the dynamical system describing the dynamics of our cosmological model 
is given by the following set of equations 


P 


, n k 

^ > -r 

Jn 


X' 

where q = such that 

where T? = + Yff. 


M{t) + X{l + q-ml) , 

s±(g - 2), 

+X2 = 1, 


(24) 

(25) 

(26) 


3 The Conditions for Eternal Inflation 

We will now derive conditions for the Bianchi Type I model considered in this paper to 
exhibit eternal inflation. In the non-eternal, slow-roll inflation scenario, one defines a 
slow-roll parameter, e as: 

H' 

Comparing with the general definition for the deceleration parameter q |44j . 

H'=-{I+ q)H, 

we see that slow-roll inflation will take place whenever 

e < 1 =;> g < 0. (29) 

However, because we are seeking conditions for eternal inflation, this condition may not 
be sufficient. As discussed in [48] and [49], for inflation to be eternal, the background 
spacetime requires frequent violations of the Null Energy Condition. One can typically 
take fluctuations of the energy-momentum tensor of the scalar field as a possible source of 
such violations. In the situation where quantum fluctuations dominate slow-roll evolution, 
the energy-momentum tensor fluctuations are large enough to create Null Energy condition 
violations. Thus, in this section, we shall derive conditions where the quantum fluctuations 
dominate the slow-roll evolution. 


(28) 
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To derive the conditions for eternal inflation, we first note that a stochastic differential 
equation has the general form 


dX = a{X,t)dt + b{X,t)dW, 


(30) 


where a{X,t) is known as the drift coefficient which effects the deterministic evolution 
of the system, while b{X, t) is known as the dispersion coefficient which introduces ran¬ 
dom/stochastic fluctuations into the system. Using Eq. (|19ll . we can write Eq. (I24p as 

dX =[{l + q-ml)X]dT + dW. (31) 


Comparing these two equations, we see that the dispersion coefficient is unity, while the 
drift coefficient is given by (l + (7 — tttq) X. Eor eternal inflation to occur, the quantum 
fluctuations which are represented by the stochastic forcing term in Eq. (|3ip must dominate 
over the deterministic slow-roll term [39]. In a stochastic differential equation of the form 
Eq. (pup , the dynamics of X will be dominated by stochastic ‘jumps’ if b{X,t) > a{X,t). 
Therefore, to have eternal inflation, we require that 

{l + q-ml)X <1, (32) 

where because of the Eriedmann equation (1261) . we have that —1 < X < 1. The restriction 
given by (I32p will further constrain (1291) as we will now see. The only way to satisfy (I32|) 
is to take _ 


-1 < A < 0 , 

0 < mo < Y 

-l + X + qX 

X 

(33) 

- 2 < g < 0 , 

0 < A < 1 , 


(34) 


where, of course, mo > 0. Therefore, as can be seen from (l33P and (f33p . the condition 
()29p is not sufficient for eternal inflation. Taking — 2 < g < 0 in combination with the 
conditions in (|33p and (|34p will ensure that eternal inflation takes place. 

One can also go a step further by writing according to Eq. (1261) X"^ = 1 — T? = 
1 — . Making this substitution in (j33|) and (I34p . we obtain the following restrictions 

on the anisotropy of the model: 


V3 


< 


a/3’ 


V3 


< s_ < 


\l 


1-S/ 

V3 


(35) 


An equivalent expression can also be obtained by interchanging S_|_ and S_ in (|35p . 

One sees that the anisotropy required for eternal inflation to take place must be con¬ 
tained in a small region within the Kasner circle, -|- S/ = 1. We have displayed the 
regions of anisotropy that allow for eternal inflation in Eigure [H based on the conditions 
(I35p . These conditions support the claim made in [7|, where inflation will only occur if there 
is not a significant amount of anisotropy to begin with. However, we believe this is the 
first proof using expansion-normalized variables based on an orthonormal frame approach. 








I- 



Figure 1: A plot of the Kasner circle = 1. For eternal inflation to occur, the 

anisotropy variables must take on values (S_|_,S_) within the shaded region. 

4 The Probability of Eternal Inflation 


We now wish to calculate the probability of eternal inflation occurring. We first note that 
Eqs. (j24p - (j25p upon using the definition in Eq. (jl9l) take the form 


dX = [{1 + q - ml) X] dr + dW, 

(36) 

cT 

1 

-H 

w 

II 

-H 

w 

(37) 

where q = 2S^ — A^, such that 


+ A^ = 1, 

(38) 

where 



Because of the symmetry of the Bianchi Type I model, from the Friedmann equation, 
we will write = 1 — A^, which implies that q = 2 — 3X^. So, the inflaton dynamics is 
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governed by a single stochastic differential equation 


dX =[{l + q- ml) X] dr + dW, 

(39) 

where 


q = 2- 3X^. 

(40) 

In fact, it is convenient to substitute q into Eq. (I3^. and write 


dX = [X (3 - ml) - 3X^] dr + dW. 

(41) 


Because of the highly nonlinear drift term in Eq. (1411) . there is no closed-form solution 
to Eq. m [50] . In fact, what is often of interest in applications of stochastic differential 
equations is the expected value of the process X, which in this case, also has no closed-form 
expression. Another interesting thing to note is that the deceleration parameter q is now 
a random variable, as can be seen from Eq. (14011 . where X{t) is given by Eq. (1411) as 


X{t) 


X{0)+ nxis) {3-ml) 

Jo 

X{0)+ Hxis) {3-ml) 

Jo 


3X{s)^]ds+ [ dW{s) 

Jo 

3A:(s)3] ds + W{T). 


(42) 


Eurther, recall from the previous section that the condition for inflation to take place 
is given by the conditions (j33p and (I34p . Erom Eq. (j40l) . we see that inflation will occur 
when 


l<x<i^^l< 


A(0) + / [X(s) (3 - ml) - 3X(s)3] ds + W(t) 


< E 


(43) 


or 


- 1 < X < 


^ -1 < 


X(0) + / [X(s) (3 - ml) - 3X(s)3] ds + W{t) 


< 



where in this latter equation, we have the condition on mo as described in (j33p . 

We are now in a position to calculate the stationary probability distribution. The sta¬ 
tionary distribution is of great interest with respect to eternal inflation, since the stationary 
distribution describes the distribution of X (r) over a long period of time. In addition, it is 
well-known that inflation should last for more than 60 e-folds [2]. We proceed as follows. 
Eollowing m, we first compute the scale measure, s(X), as 


s(X) 


exp 


exp 



- rnl) - 3y^dy 



Mi, nM .. Mo ^ 

2 2 4 J 


(45) 
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The speed measure, m{X), is then computed as 


m{X) = 


s(X) 


3X2 ^2^2 3J^4 


= exp 

The stationary density, 7r(X), is then given by 

m{X) 


^^0 3 Xq 


( 46 ) 


7r(X) = 


m(X)dX 

^_lx2(2m2+3X2-6) 


(47) 


The motivation behind choosing the limits of integration as equal to plus or minus unity 
for the stationary density calculation in Eq. (|171l is the Friedmann equation, Eq. (l38t) . 
which tells ns that X = 0 when S2 = 1, and X = ±1 when Tp‘ = 0. 

It is perhaps of some interest to calculate this probability of eternal inflation occurring 
using some values for mg. Following m, m, and [53], one can take H ps 10^2 QeV 
(10“^ in natural units), and m k, 10^2 _ GeV (10“^ — 10“® in natnral units). We 
shall therefore calculate the probability of eternal inflation occurring for two values of ttiq. 
Recalling the definition of mo from Eq. (I16p . we see that uiq = for m = 10^2 QeV, 

and niQ = ^, for m = 10^^ GeV. 

v3 

First, for mn = -)=, we calculate 
v3 


P(X) = 




= 0.299. 


0.253641e"^^"(^^""^)dX + [ 0.253641e"5^"(^^""^)dX 


(48) 


Second, for mn = we calculate 
v3 

P(X) = 




hlll09e"^^"(^^"+^)dX 


= 4.87 X 10 


-11 


(49) 


The BIGEP2 experiment recently determined that H ~ 10^^ GeV (10 ^ in natural 
units) during inflation |5l]. Using this value, we now try calculating the eternal inflation 
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probability with tuq = and mo = We obtain that for mo = 


P{X) = 


r \/2/3 lv2^ov2 89999 ^ 1 v2 ^ v2 89999 ^ 

0.21677e"2'^ -JmoJdX + / 0.21677e"2'^ -JsmJdX 


J-i 

= 0.3367. 


'VW 


(50) 


While for mo = we see that 

P{X) = 




0.217118e"^'^"(^^""is)dX+ / 0.217118e"5^"(^^""in)dX 


J-i 

= 0.3363. 




(51) 


We thus see that the highest probability of eternal inflation occurring, about 33.67% 
occurs for mo = l/(100\/3), which corresponds to m = 10“^, and H = 10“^. The second- 
highest probability of eternal inflation occurring was 33.63% for mo = l/(10\/3), which 
corresponds to m = 10“® and H = 10“^. We also see that the probability of eternal 
inflation is infinitesimally small when mo = 10/\/3, which corresponds to m = 10“® and 
H = 10-'^. 


5 Conclusions 

In this paper, we studied the phenomenon of stochastic eternal inflation for a chaotic infla¬ 
tion potential in a Bianchi Type I spacetime background. After deriving the appropriate 
stochastic Klein-Gordon equation, we derived the conditions for eternal inflation to occur. 
We showed that for eternal inflation to occur, the amount of anisotropy had to be very 
small. In fact, we showed that eternal inflation would only take place if the shear anisotropy 
variables took on values within a small region of the interior of the Kasner circle. We then 
calculated the probability of eternal inflation occurring based on techniques from stochastic 
calculus. 

We found that the highest probability of eternal inflation occurring, about 33.67% 
occured for mo = l/(100\/3), which corresponds to m = 10 and H = 10 The 
second-highest probability of eternal inflation occurring was 33.63% for mo = l/(10\/3), 
which corresponds to m = 10“® and H = 10“^. We saw that the probability of eternal 
inflation was infinitesimally small when mo = lO/y/O, which corresponds to m = 10“® and 
H = 10-^. 
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